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We argue that the large N strong/weak phase transition is a generic phenomenon in a 
finite temperature supersymmetric Yang-Mills theory of maximal supersymmetry. M = 4, 
D = 4 SYM is the canonical example, where we also argue that the large N Hawking- 
Page phase transition disappears for a sufficiently small coupling. The Hawking-Page 
transition temperature is lowered by the first a' correction. Physically, the strong/weak 
phase transition is identified with the correspondence point of Horowitz and Polchinski. 
We also try to construct toy models to demonstrate the large N phase transitions, with 
limited success. 
Oct. 1998 



1. Introduction 

Supersymmetric Yang-Mills theories (SYM) exhibit rich dynamics. The study of ma- 
trix theory and Maldacena's conjecture |||| reveals much of the unusual physics of 
SYM with maximal supersymmetry in the large N limit. Small energy gaps, different en- 
ergy scales are only a few examples of many intriguing properties of these theories. While it 
is surely a great hint that both matrix theory and AdS/CFT correspondence encode grav- 
itational physics in the large N limit, the major technical obstacle to our understanding of 
these models is to truly resolve the large N problem. 

Another surprise the AdS/CFT correspondence brings to us is the possibility of solving 
the confinement problem of QCD in the large N limit. The hope might be hindered by 
a large N strong/weak phase transition [|J. According to ||, D4-brane and D3-brane 
theories when compactified on a Euclidean circle lead to QCD and 3D pure Yang-Mills 
theory, provided the radius of this circle is much smaller than the intrinsic scale of the lower 
dimensional theory. The latter can be viewed as the low energy sector of the original theory 
at a finite temperature. The strong/weak phase transition of the free energy reflects the 
same physics of this low energy sector in the D = 4 case, since its occurrence is independent 
of the temperature. One can not completely ignore the underlying 4 dimensional theory 
at the transition point, since the 3D YM scale A3 = AT ~ T, the KK modes have the same 
energy scale as the 3D QCD string tension. While in the D = 5 case, although the phase 
transition is in terms of temperature, so long when the M-circle scale is small enough, 
the transition temperature can be very high, again the phase transition will be reflected 
by physics of the low energy sector, namely QCD. The existence of the phase transition 
implies that to study the weak coupling regime, one must study the weak coupling regime 
directly. The weak coupling physics can not be obtained by the analytic continuation of 
the strong coupling physics. 

The study of |4[] was conducted in the infinite volume limit, and is limited to D = 4 
case. We shall generalize it to the finite volume case, and point out that the strong/weak 
phase transition also exists on S 3 . The phenomenon is also a general feature of SYM of 
maximal supersymmetry in various dimensions. The phase transition point is identified 
with the correspondence point of || . This makes the definition of the correspondence point 
precise. For the p + 1 dimensional SYM, if p 7^ 3, we predict the phase transition temper- 
ature T c ~ A 1 ^ 3- ^. Our study of phase transition is rather qualitative. To understand 
dynamic details, one need to study both the weak coupling side and the strong coupling 
side to demanding precision. There is no tool available to enable us to do so. 
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We shall point out the connection between the strong/weak phase transition and the 
correspondence principle in the next section. In the subsequent section, we argue that the 
large N first order Hawking-Page phase transition disappears in the weak-coupling regime 
(D = 4), thus there is a phase transition of phase transition, which is just the strong/weak 
phase transition. The correction to the Hawking-Page transition temperature due to the 
first a' correction is found in sect. 4, where we determine the critical coupling A c using this 
result. In the final section we design some toy models to explain large N phase transitions. 
We show it is quite hard to come up with a good model to demonstrate a first order phase 
transition. We point out possible missing physics in these models. 



2. Strong/weak phase transition and correspondence principle 

It was argued in [Q that the free energy of M = 4, D = 4 SYM at a finite temperature 
is not a smooth function of the coupling constant A = (7y M iV. To determine the nature 
of this phase transition, such as its order, is beyond the power of the currently available 
tools. 

Denote the critical value of A by A c . For a A < A c , one can trust the SYM pertur- 
bative expansion. In this regime, one does not expect the origin A = be an essential 
singularity, given the fact that the underlying zero-temperature theory is finite, and the 
instanton effects are suppressed in the large N limit 0. On the other hand, for a A > A c , 
one need to use the whole string theory on manifold AdS§ x S 5 . This reminds us the 
correspondence principle of Horowitz and Polchinski || . They postulate that there exists 
a correspondence point at which the maximal curvature invariant at the horizon becomes 
of the string scale. For a smaller curvature, the system is in its black hole phase, and the 
semi-classical geometry is reliable in string theory. One can study the closed string theory 
in this background. For a larger curvature, the geometry loses its meaning in string theory, 
now the D-brane perturbation theory takes over. Our specification of a phase transition 
defines precisely what the correspondence point is. 

Here we must note that the following argument is not meant to prove the existence 
of a phase transition, rather, we intend to argue that if there exists a phase transition, it 
is natural to identify it with the correspondence point. This identification has proven to 
be useful in recent work [11,12]. 

The connection between the strong/weak phase transition and the correspondence 
principle looks puzzling at the first sight. For a standard Schwarzschild black hole, the 
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curvature invariants would be proportional to 1/Vq, where rg is the horizon size. This 
quantity depends on the Hawking temperature. We know that in a conformally invariant 
theory with an infinite volume, there can be no phase transition in terms of temperature. 
The fact that we are applying the correspondence principle to a AdS black hole saves the 
day. We will show, that all the curvature invariants in such a background are independent 
of 7"o, so long as tq > 0. Thus the maximal curvature is determined by the AdS size, 
R = (2A) 1/4 / s . Now R ~ l s if A ~ 0(1). This is consistent with the fact that A c ~ 0(1). 
The metric of the AdS 5 black hole, with a boundary R x S , is @ 

*> 2 = - iw 2 +!>?] + - ir 1 ^- (2.1) 

i=l 

With this metric, it is not straightforward to see that all the curvature invariants are 
independent of Tq. To do so, we follow M, to rewrite the metric in a different coordinates 
system. Introduce 

r 4 = r 4 (l + 6 2 ), t=^V, (2-2) 

ir 

the metric becomes 

J^ + vTT^i:^, (2.3) 

where we also rescaled Xi. The metric in the parentheses is universal, thus its curvature 
invariants are universal too. We conclude that the curvature invariants of the AdS black 
holes are independent of ro, and determined only by the AdS size R. 

The universal metric (|2.3|) is applicable only when ro > 0. The AdS black hole solution 




is certainly different from that of AdS. The coordinates transformation ( |2.2| ) is singular in 
the ro = limit. In other words, the metric (|2.3| ) describes the very near-horizon region 
in the ro —>■ limit. 

The observation on the connection between the large N strong/ weak phase transition 
and the correspondence point generalizes to other SYM with the maximal supersymmetry 
in various dimensions. For D 7^ 4, the YM coupling constant is dimensionful, and the 
construction of a dimensionless constant involves both the coupling constant and the energy 
density. Thus the strong/ weak phase transition becomes of a finite density or a finite 
temperature phase transition. Let A = g\ M Nd v , where d p is a numeric factor depending 
only on p. The near-horizon metric of a nonextremal Dp-brane is || 



ds = a 1 



I fdt 2 + - ^_ J2 dx * + f~ ldJj2 + ^U^/ 2 dnl_ p \ , (2.4) 



where 

jt(7-p)/2 / U 7 ~ P \ 

^ — i 1 -^)- (2 ' 5 » 

The dilaton field is given by 

= (2n) 2 -?g 2 YM . (2.6) 

For p < 3, the string coupling becomes stronger when Uq is smaller. For p > 3, it becomes 
weaker when Uq is smaller. 
As in the p = 3 case, use 

U 7 -? = U 7 Q - p (l + b 2 ), t=^—U^^f s ^ (2.7) 
where A e // = A£/q -3 , the metric is put into 

2 Aa' J\ e ff ( 2 3 P -17 „ 6 2 rf^ 2 (7-p) 2 ^ ,2: p ~ 3 \ 

ds = (7 _p)2 (J 1 + 6 ) 14 " 2P rf6 + 7f=p + ^"T^ C 1 + 6 ) 14 - 2P J 



We see that the curvature is controlled by X e ff- When this quantity becomes order 1, the 
maximal curvature becomes of the string scale. Since the Hawking temperature is 

T = (A e// )" 1 / 2 t/o(^), (2.9) 

the temperature can be identified with Uq at the correspondence point. This in turn 
implies that the effective coupling at the correspondence point is A e // = AT P_3 ~ 1, and 
the phase transition occurs at 

T C ~A^. (2.10) 

The dilation field at the horizon is 

e^^(Xeff) 7 -^. (2.11) 

it is of order 1/N at the phase transition point. The large N approximation becomes quite 
good, and on the string side one is studying free strings propagating in the black hole 
background. 



When the above analysis applied to p = 4 case, we predict that the phase transition 
occurs at T c ~ A -1 . For higher temperatures, A e // is larger, and we are in the supergravity 
phase. For lower temperatures, we need to apply the perturbation theory of SYM or the 
underlying M5-brane theory. To study QCD following ||, temperature T must be adjusted 
over the induced QCD scale T ^> Aqcd- At the QCD scale, the four dimensional QCD 
coupling A4 = AT receives quantum corrections, but it is still of order 1. Thus, the situation 
is close to the transition point. We believe that it is located at the small coupling side. 
One of hints is the observation that the Regge trajectory in the strong coupling regime 
displays very different pattern as expected of QCD 0. 

The 5 dimensional Yang-Mills coupling is given by g\ M = Ru, where Ru is the radius 
of the M-circle around which the M5-branes are wrapped to obtain D4-branes. Thus the 
phase transition occurs at (3/Rn ~ N. 

Another interesting case is p = 1. The transition temperature is given by T c ~ A 1 / 2 . 
For T ^> T c , the theory is strongly coupled. If we switch the role of time and space, the 
theory becomes the 2D SYM living on a circle of circumference f3. This can be related 
to matrix string theory ||10|| . In the strong coupling regime, the theory is effectively a 
conformal orbifold model, and the physics is dominated by long-strings. For 1/T small 
enough, one can ignore fermions, since fermions are anti-periodic. The theory is effectively 
bosonic. When T is lower than T c , the theory is nonabelian. The phase transition separates 
the abelian orbifold phase of the bosonic sector and the nonabelian phase of the NS sector. 

We turn to AdS$ with boundary S 3 x S 1 in the next two sections. For discussions on 



the phase diagrams with other topologies, see ['I !.] and [12 



3. Disappearance of Hawking-Page phase transition 

Hawking and Page found some time ago [|13|], that there are two manifolds for a 
given temperature on AdS space, provided the temperature is not too low. This was re- 
interpreted by Witten as a large N phase transition in the Af = 4 D = 4 SYM |5| . This is 
a finite volume effect, and occurs only when the SYM lives on S 3 x 5" 1 or S 2 x R x S . 

The AdS black hole metric with boundary S 3 x S 1 is 

ds 2 = (1 + ^ " ^ + (1 + ^ - ;!)" V + r^dnl (3.1) 

For a given Hawking temperature, there are two black holes, one small, one large. The 
smaller one can not be regarded as a thermal state in the boundary Yang-Mills theory, 
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since its entropy is not an extensive quantity. The horizon radii of these two black holes 
are 

Use the standard formula for the Bekenstein-Hawking entropy A/4G§, the entropy 
reads 

(3.3, 

where we used formulas G5 = Gio/tt 3 R 5 , and Giq = 2 3 ir 6 g 2 (a') 4: . In the large volume or 
high temperature limit, S- is not an extensive quantity. Compare S+ with the one in the 
infinite volume limit |fT3| , we find that V3 = 2n; 2 R 3 , namely the radius of S 3 is also R. 
Note that the precise definition of the radius of S 3 depends on the choice of time, or that 
only TR is convention independent. 

There is a finite size correction to the entropy formula. And formula ( |3.3| ) displays 
a branch cut at ivTR = y/2. There is a minimal Hawking temperature. However, the 
Hawking-Page transition temperature is tcT c R = 3/2. This value is obtained by demanding 
the subtracted action vanishing, thus = R 2 . The critical temperature is higher than 
the minimal temperature, thus there is no particular physics associated with the minimal 
temperature. 

The free energy can be obtained by using the relation S = (3 2 dpF. 

One can trust the above picture only when classical supergravity is valid, thus the 
radius of S 3 is much larger than the string scale. In terms of SYM, this requires A ^> 1. 
We now argue that for a small A, the first order phase transition disappears. Denote the 
phase transition temperature by T C (X), it is a function of A. In the large A limit, T C R 
tends to a constant, 3/(2%). If the phase transition disappears for small A, there must 
be a critical A c , at which T c = 0. This phase transition of phase transition must be the 
strong/weak phase transition. 

Apparently, T c = at A = 0. The theory is free, and the free energy can be computed 
by calculating determinants of free fields. The number of free fields is of order 0(N 2 ), so 
the free energy scales as F ~ A 2 , and we always have the high temperature phase. For 
a small A, to include effects of interaction, one employs the loop expansion in A, and in 
the large N limit, it is sufficient to count the planar diagrams. At each loop order A n , 
the number of planar diagrams grows as c n [jl5|jy| . There is no UV divergence, since the 
underlying zero temperature theory is finite. There is also no IR divergence, since both 
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the space and the time have a finite size. Thus, the systematic loop expansion of the free 
energy assumes the following form 

F = -N 2 T A R 3 A n /n(Ti?), (3.4) 

n=0 

where f n depends only on the dimensionless combination TR. We do not expect effects 
such as UV or IR renormalons, the above series ought to be convergent for a small enough 
A, and the series typically has a finite convergent radius. 

By a similar argument, each perturbative series at a given genus also has a finite 
convergent radius. Now it is possible that a singularity is located on the positive axis of A, 
denote this by A c . Beyond this point, one can no longer trust the gauge theory perturbative 
expansion, and one must invoke other techniques such as the closed string theory on the 
AdS space to compute the free energy, for A < A c , the free energy always scales as N 2 , and 
there is no Hawking-Page phase transition. The critical A c is determined by the asymptotic 
behavior of f n (TR). It is still possible that A c is a function of TR. Consider the limit 
TR — > oo. In such a limit, as argued in jjjj, there is a finite critical A c . When one lowers 
TR, it is possible that A c is also lowered, but cannot be pushed all the way to zero. For 
instance, assume that the convergent radius is proportional to T a with a positive a, then 
the asymptotic form of f n ~ {l/TR) na . Analyzing the Feynman diagrams, such a factor 
can come only from the propagators of the zero modes. Now on S 3 , there is no zero mode 
at all. Put in another way, there is a mass gap in the SYM on S 3 which is proportional to 
the smallest scaling dimension. Still denote the minimum of X C (TR) by A c . We conclude 
that there is no large N Hawking-Page phase transition below A c . 

The existence of Hawking-Page phase transition indicates that there exists a finite A c , 
above which the perturbative SYM ceases to be valid, and the AdS picture takes over. We 
shall show in the next section that the Hawking-Page temperature is lowered by the first 
a' correction. 



4. ol correction to the thermodynamic quantities 



The correction to thermodynamics by the first a' correction was done in | 17|j , in the 
infinite volume limit, where it was found that the correction to the entropy is positive. 
There are two ways to compute this correction. The simpler way is to substitute the Rie- 
mann curvature of the metric ( |3.1| ) directly to the term (o/) 3 /? 4 to compute the subtracted 



action. Another is to take the corrected metric, dilaton into account. The two calculations 



give the same result. The full answer of the corrected metric was found in [[18|| , again in 
the infinite volume limit. 

We will adopt the simpler method, hoping that the corrected metric will not change 
the final results. The full form of the (a') 3 i? 4 is 



AI 



^/^C(3)(a') 3 ^ + ---, (4-1) 



16nG 

where the dots denote terms depending on other fields, and 

^pmnq^h ^rsk + 2 ^pqmn^ h ^ rsk - \^-^) 

Introducing angular variables x, <fi on S 3 , after some calculations, we obtain those 
nonvanishing components of the Weyl tensor 



° cd 


= 3Xe ab e cdl a,b = t, r, 




r~iai 
° bj 


= -X8t8% i,j = X ,e,(t>, 


(4.3) 


° kl 


= X(5lSi-5i8i), 





where X = r^/r 41 . The difference between our result and that of [17| is that X is indepen- 
dent of R 2 in our case. 

It is straightforward to substitute the above result into (|4.1|) , the result is 



15ttC(3) (a'frlfi 
~64 G 5 r^ 



AI=- 2 y ^ J 2 ■ (4-4) 



Using the formula AS = P^dpAF we obtain 
AS = 157r2C(3) iV 2 (2A)- 3 / 2 y 3 T 3 [- - -(1 - _?_)i/2]3[3 _ (i _ 2 



-1/21 



32 v ' ° L 2 2 V (ttTR) 2 ' 1 1 K (tvTR) : 



(4.5) 



In the large volume limit we recover the result of [|17| . Note that the last factor approaches 
minus infinity when T approaches its minimal value. At the uncorrected transition tem- 
perature T c , it becomes zero. 

Demanding the subtracted action be zero, we calculate the corrected Hawking-Page 
temperature, up to the leading order (en') 3 

TiZTc = \ - ^C(3)(«') 3 AR 6 = I - Y C(3)(2A)" 3 / 2 . (4.6) 



As expected, the corrected temperature is lower. The correction to the entropy (O) at 
this point is weighted by A -3 . The last factor will never be too negative to render the total 
entropy become negative. 

Although the critical temperature is lowered, we can not conclude that T c will drop 
to zero at a finite A c . If we are optimistic and assume the higher a' terms can be ignored 
at A c , Eq. (|4.6|) can be used to determine A c : 

A c = ^(5C(3)) 2 / 3 + .... (4.7) 

This value is greater than 1. 

As in sect. 2, the connection between the strong/weak phase transition point and the 
correspondence point is still valid. For the larger black hole, since r_|_ > R/y/2, the maximal 
curvature is always controlled by R. The smaller black hole is more like a Schwarzschild 
black hole. Its specific heat is negative, and can be described at best as a meta-stable 
coherent state in the SYM theory. The horizon can be arbitrarily small, and the maximal 
curvature at the horizon, unlike the case discussed in sect. 2, can be arbitrarily large. 



5. Some toy models 

We design some artificial toy models in this section to explain the occurrence of the 
large N phase transitions. To our knowledge, the exactly solvable models such as those 



studied in and |2^] exhibit different, although similar, phase transitions to that of 
Hawking and Page. Those models are all effectively one-matrix models, and the phase 
transition is of the third order. The Hawking-Page transition is a first order one, since the 
free energy scales differently in N in the low and high temperature phases. 

Consider a conformal field theory defined on spacetime S X R. Let H denote the 
Hamiltonian on S d . The free energy is given by 

F = -ln(]Te-^). (5.1) 

The metric on S x R reads 

ds 2 = dr 2 + R 2 dtt 2 d , (5.2) 

where R is the radius of S d . This metric is conformal to a flat metric on R d+1 . Use the 
new coordinate r = exp(r/R), the metric is mapped to 

ds 2 = % (dr 2 + r 2 dn 2 ) , (5.3) 
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the new metric in the parentheses is the flat one on R d+1 . Now the time translation 
r — > r + a gets translated into rescaling of r: r —>■ exp(a/R)r. Thus the Hamiltonian on 
S d is identified with D/R, where D is the dilatation operator on R d+1 acting at the origin. 
For simplicity, we still use f3 to denote the combination f3/R. The task of computing the 
free energy now becomes of counting scaling dimensions: 

F = -In QTe"^). (5.4) 

A large N quantum field theory is very complicated, so we will construct simplest 
possible large N models, one matrix models or multiple-matrix models. Given a Hermitian 
one matrix (p, there are only N independent single trace operators tr^ 71 , n = 1, . . . N. A 
"multi-particle" state, similar to those considered in [pQ| , is given by the product of a 
string of such single trace operators, 

tr(P ni ...tr(P ni . (5.5) 

Assign a scaling dimension A(n, A) to operator tr</> n , where we introduced a coupling 
constant A. The scaling dimension of a multi-particle operator is generally complicated. 
To simplify the situation, we assume the scaling dimension of the general operator be given 
by A(?ii, A). The partition function is then 



N 



n=l 

where q = exp(— /?). The free energy 



'N 

n=l 



= J> 1-^ . (5.6) 



The first model is specified by A(n, A) = n. All operators mimic the chiral primary 
operators. In such a case, we expect no phase transition when the temperature is varied. 



This model looks the same as the 2D QCD string on a torus |19[]. The free energy is 
calculated by expanding the logarithmic 



00 k 00 k 



k ' f^> k 

fc=i fc=i 



There is no 1/N corrections. The second sum is due to the instanton contribution. Each 
term can be interpreted as the k-instanton contribution. 
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Our second model is given by A(n, A) = a(A) Inn + b(X) with b > 0. If a(A) 7^ 0, for 
large n, the scaling dimension receives large quantum correction. These operators mimic 
operators corresponding to stringy states. 6(A) > is to ensure the positivity of the first 
scaling dimension. 

Again, expanding the logarithmic in (|5.6|), we have 



00 „-kj3b N 
fc=l n=l 

Apparently, the behavior of the finite sum depends crucially on the value of (3a. The 
dominant contribution to Fn comes from k = 1. 

2^ n In N, (3a = 1 1 J 

n=l 

The model exhibits kind of phase transition. When (3a < 1, or T > a, the free energy 
scales as a positive power of N; it tends to a constant when T < a, and grows as lnA^ at 
T = a. So T c = a is the critical temperature. Unfortunately, there is no standard large N 
expansion. To see the nonanalyticity in N, we write 

00 -kf3b 

F " = ~Y1 —j^-[Wa) - ((k(3a, N + 1)] = /(l) - /(TV + 1). (5.10) 
fc=i 

the first term is independent of N, and is regular provided (3a 7^ 1/k. To examine the 
second term, use the Hermite formula 



C(k(3a, N + 1) =l(N + l)- k(3a + io l A N + 
2 k(3a — 1 



(5.11) 



+ 



/■oo 

2 / dy((N + l) 2 + y 2 )-W 2 sm(k(3a6)(e 27Ty - l)" 1 , 

where 6* = arctan(y/(iV + 1)). One finds 

/(TV + 1) =f crit (N + 1) + (N + 1/2) ln[l - (N + l)-^e"^] 

/■oo 

- / dy(e 27Ty - l)" 1 ln{[l - (y/(N + l) 2 + y 2 e - ia0+b )- ] (5.12) 

[1 - ( v / (iV+l) 2 + 2/V ae+& )- /3 ]- 1 }, 

where the first term and the second term control the large N behavior. The first term is 
given by 

OO n 

f crit (N + 1) = J2 k^\ e ~ kpb{N + l)1 ~ kPa - (5 - 13) 
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We are convinced that it is impossible to construct a one-matrix model to exhibit 
both a first order phase transition and the usual large N expansion. This leads us to 
examine multiple-matrix models. As soon as there are more than one matrix, the number 
of single-particle states at a given level starts to proliferate. For instance, consider two 
matrices. A general single-particle state is 

tr(^r^--^r^)- (5-14) 

Define the level n = XX m « + n i)- At a given level, the number of the single-particle states 
grows at least as the partition number P(n), since there is no longer permutation symmetry 
among m; and n*, only the cyclic symmetry is retained. This growth of the single-particle 
states is the generic property of a string theory. 

For finite N, one must implement constraints on multiple-particle states. These con- 
straints are known as the Mandelstam constraints, and as the stringy exclusion principle 



in the modern context [20|. It is difficult to implement these constraints efficiently. We 
will introduce a cut-off in the same fashion as in (|5.6|) . The class of our multiple-matrix 
models is specified by 

N 

Z N = Y[ (1 - g A(n,A))-p(rO (5.15) 
n=l 

where we assumed that all the single-particle states at a given level n have the same scaling 
dimension A(n, A), and the multiplicity is given by p(n). Needless to say, these models are 
over-simplified multiple-matrix models. 

We show that some physics is missing in these toy models, so that it is impossible 
to come up with a first order phase transition. First, we take A(n, A) = a(X)^/n + b(X). 
The factor yfn is introduced according to 0. As we argued before, the growth of p(n) 
is at least close to the partition number P(n), thus p(n) ~ exp(cy/n). This introduces a 
problem, namely there exists a limiting temperature, the Hagedorn temperature. Above 
the Hagedorn temperature, the free energy will increase as an exponential of y/N. This is 
not allowed. So we shall assume that due to Mandelstam constraints, the growth of p(n) 
is weaker than exp(c\/n). 



Take the logarithmic of (|5.15|) , the sum can be approximated by an integral, in the 
large N limit. The upper bound of the integral can be pushed to infinity, since we assumed 
that there is no limiting temperature. Thus, the free energy is given by, in the large N 
limit 

F = J (5.16) 
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where we assumed that p(x, N) in general is a function of N too. If there were a critical 
/3 C at which a first order phase transition happens, the first derivative of F would diverge. 
Take the first derivative 



We now require that the above integral diverge for a (3 C . The divergence should not come 
from x = oo. If it comes from x = xq, then xq must be a singularity of p(x, N), since the 
other factor is regular everywhere. In this case, the integral will diverge for all /?, not just 
for a single (3 C . 

We conclude that our simple-minded toy models do not capture the physics of Af = 4, 
D = 4 SYM at a finite temperature. Many things might have gone wrong in our models. It 
is possible that solution to the Mandelstam constraints can not be mimicked by our ansatz. 
It is also possible that the assumption that the scaling dimension of a multiple-particle 
operator is just A(njA) is badly wrong. It remains a challenging problem to construct 
a more realistic model to demonstrate both the first order Hawking-Page phase transition, 
as well as the large N strong/weak coupling phase transition. 
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